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1. Introduction 

Nowadays the question 'what is complexity?' is a challenge to be answered. This question 
is triggering a great quantity of works in the frontier of physics, biology, mathematics and 
computer scienc e. Even more when this century has been told to be the century of complexity 
(Hawking, 2000). Although there seems to be no urgency to answer the above question, many 
different prop osals that have been developed to this respect can be found in the literature 
<Perakhl.l200l . In this context, several articles concerning statistical complexity and stochastic 
processes are collected in this chapter. 

Complex patterns generated by the time evolution of a one-dimensional digitalized coupled 
map lattice are quantitatively analyzed in Section [2] A method for discerning complexity 
among the different patterns is implemented. The quantitative results indicate two zones 
in parameter space where the dynamics shows the most complex patterns. These zones are 
located on the two edges of an absorbent region where the system displays spatio-temporal 
intermittency. 

The synchronization of two stochastically coupled one-dimensional cellular automata (CA) is 
analyzed in Section [3] It is shown that the transition to synchronization is characterized by 
a dramatic increase of the statistical complexity of the patterns generated by the difference 
automaton. This singular behavior is verified to be present in several CA rules displaying 
complex behavior. 

In Sections H] and |5l we are concerned with the stability analysis of patterns in extended sys- 
tems. In general, it has been revealed to be a difficult task. The many nonlinearly interacting 
degrees of freedom can destabilize the system by adding small perturbations to some of them. 
The impossibility to control all those degrees of freedom finally drives the dynamics toward 
a complex spatio-temporal evolution. Hence, it is of a great interest to develop techniques 
able to compel the dynamics toward a particular kind of structure. The application of such 
techniques forces the system to approach the stable manifold of the required pattern, and then 
the dynamics finally decays to that target pattern. 

Synchronization strategies in extended systems can be useful in order to achieve such goal. In 
Section |U we implement stochastic synchronization between the present configurations of a 



cellular automata and its precedent ones in order to search for constant patterns. In Section[5] 
this type of synchronization is specifically used to find symmetrical patterns in the evolution 
of a single automaton. 

2. Complexity in Two-Dimensional Patterns Generated by Coupled Map Lattices 

It should be kept in mind that in ancient epochs, time, space, mass, velocity, charge, color, etc. 
were only perceptions. In the process they are becoming concepts, different tools and instru- 
ments are invented for quantifying the perceptions. Finally, only with numbers the scientific 
laws emerge. In this sense, if by complexity it is to be understood that property present in all 
systems attached under the epigraph of 'complex systems', this property should be reasonably 
quantified by the different measures that were proposed in the last years. This kind of indi- 
cators is found in those fields where the co ncept of information is crucial. T hus, the effective 
meas ure of complexity iGrassbergerl 1 19861) and the thermodynamical depth jLlovd & PagelsL 
19881) come from p hysics and other attempts such as algorithmi c complexity iKolmogorovt 
19651; IChaitinl. 1 1966|) . Lempel-Ziv complexity jLempel & Zivl.ll976l) and e-machine complexity 
JCruthfieldl. 1 1989) arise from the field of computational sciences. 

In particular, taking into account the statistical properties of a system, an indicator called 
the LMC (LopezRuiz-M ancini-Calbet) complexity has been introduced jLopez-Ruid. 1 19941: 
lLopez-Ruiz et aUll995h . This magnitude identifies the entropy or information stored in a sys- 
tem and its disequilibrium i.e., the distance from its actual state to the probability distribution 
of equilibrium, as the two basic ingredients for calculating its complexity. If H denotes the 
information stored in the system and D is its disequilibrium, the LMC complexity C is given by 
the formula: 

C(p) = H(p)-D(p) = 

= -k(Ll 1 P,lo S p l )-(Ll l {p t -i) 2 ) (1) 

where p = {p,}, with p,- > and i = ■ ■ ,N, represents the distribution of the N accessible 
states to the system, and k is a constant taken as 1 / log N. 

As well as the Euclidean distance D is present in the original LMC complexity, other kinds of 
disequilibrium measures have been proposed i n order to remedy some statis tical character- 
istics considered troublesome for s ome authors iFeldman & Crutchfieldl, ll998). In particular, 
some attention has been focused |Lm1. Il99ll: iMartin et all l2003h on the Jensen-Shannon di- 
vergence Djs as a measure for evaluating the distance between two different distributions 
(Pi'Pi)- This distance reads: 

Djs(pi,p2) = H(nipi + n 2 p 2 ) - 7tiH(pi) - n 2 H(p 2 ), (2) 

with tl-[,K2 the weights of the two probability distributions {p\,pz) verifying TC\,Ti2 > and 
K\ + 7t2 = 1. The ensuing statistical complexity 

C JS = H ■ D IS (3) 

becomes in tensive and also keep s the property of distinguishing among distinct degrees of 
periodicity l lLamberti et alj,l2004l) . Here, we consider p2 the equiprobability distribution and 

TTl = TZ 2 = 0.5. 

As it can be straightforwardly seen, all these LMC-like complexities vanish both for com- 
pletely ordered and for completely random systems as it is required for the correct asymptotic 
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Fig. 1. f> versus p. The /3-statistics (or BW density) for each p is the rate between the number of 
black and white cells depicted by the system in the two-dimensional representation of its after- 
transient time evolution. (Computations have been performed with Ap = 0.005 for a lattice of 
10000 sites after a transient of 5000 iterations and a running of other 2000 iterations). 



properties of a such well-behaved measure. Recently, they have been succe ssfully used to dis 
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cern situations regarded as complex in discrete systems out of equilibrium iShiner et all 
ICalbet & Lopez-Ruizll200ll;lYu & Chenl.l200d:lRosso et alJ,l2003l;l2005l;lLovallo et al.Ll2005 
As an example, the local transition to chaos via intermittency jPomeau & Mannevillelll98C 
the logistic map, x„ +1 = Ax„(l — x n ) presents a sharp transition when C is plotted versus the 
parameter A in the region around the instability for A ~ At = 3.8284. When A < Af the system 
approaches the laminar regime and the bursts become more unpredictable. The complexity 
increases. When the point A = Af is reached a drop to zero occurs for the magnitude C. The 
system is now periodic and it has lost its com plexity. The dynamical b ehavior of the system is 
finally well reflected in the magnitude C (see dLopez-Ruiz et aljjl995h . 

When a one-dimensional array of such maps is put together a more complex behavior can be 
obtained depending on the coupl ing among the units . Ergo the phenomenon called spatio- 
temporal intermittency can emerge dChate & MannevilleLll987l;lHoulriklll990l;lRolf et alj,ll998r) . 
This dynamical regime corresponds with a situation where each unit is weakly oscillating 
around a laminar state that is aperiodically and strongly perturbed for a traveling burst. In this 
case, the plot of the one-dimensional lattice evolving in time gives rise to complex patterns on 
the plane. If the coupling among u nits is modified the system can settle down in an absor bing 
phase where its dynamics is trivial {Argentina & CoulletLll99^:IZimmermann e t al, 2000) and 
then homogeneous patterns are obtained. T herefore an abrupt transition to sp atio-temporal 
intermittency can be depicted by the system iPomeaul,ll986l:lMenon et ID. l2003h when modi- 
fying the coupling parameter. 
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Fig. 2. Digitalized plot of the one-dimensional coupled map lattice (axe OX) evolving in time 
(axe OY) according to Eq. @}: if x " > 0-5 the- (z',n)-cell is put in white color and if x" < 0.5 the 
(i,n)-cell is put in black color. The discrete time n is reset to zero after the transitory. (Lattices 
of 300 x 300 sites, i.e., < i < 300 and OK 300). 



In this section, we are conc erned with measuring C and C rs in a such transition for a coupled 
map lattice of logistic type jSanchez & Lopez-Ruizll2005-al) . Our system will be a line of sites, 
i = 1,...,L, with periodic boundary conditions. In each site z' a local variable x" evolves in 
time (n) according to a discrete logistic equation. The interaction with the nearest neighbors 
takes place via a multiplicative coupling: 



= (4-3pX?K(l-x»), (4) 



where p is the parameter of the system measuring the strength of the coupling (0 < p < 1). 
The variable X" is the digitalized local mean field, 



X" = nint 



£ (*?+!+ *?- 



(5) 



with nint(.) the integer function rounding its argument to the nearest integer. Hence X" = 
orl. 

There is a biological mo tivation behind this kind of systems iLopez-Ruiz & Fournier-Prunarej. 
120041: iLopez-Ruij. 12005^) . It could represent a colony of interacting competitive individuals. They 
evolve randomly when they are independent (p = 0). If some competitive interaction (p > 0) 
among them takes place the local dynamics loses its erratic component and becomes chaotic 




Fig. 3. (•) C versus p. Observe the peaks of the LMC complexity located just on the borders of 
the absorbent region 0.258 < p < 0.335, where /3 = ( x ). (Computations have been performed 
with Ap = 0.005 for a lattice of 10000 sites after a transient of 5000 iterations and a running of 
other 2000 iterations). 



or periodic in time depending on how populated the vicinity is. Hence, for bigger X" more 
populated is the neighborhood of the individual i and more constrained is its free action. At 
a first sight, it would seem that some particular values of p could stabilize the system. In fact, 
this is the case. Let us choose a number of individuals for the colony (L = 500 for instance), let 
us initialize it randomly in the range < X; < 1 and let it evolve until the asymptotic regime 
is attained. Then the black/white statistics of the system is performed. That is, the state of 
the variable x,- is compared with the critical level 0.5 for i = 1,...,L: if x, > 0.5 the site i is 
considered white (high density cell) and a counter N u , is increased by one, or if x,- < 0.5 the 
site i is considered black (low density cell) and a counter Nj, is increased by one. This process 
is executed in the stationary regime for a set of iterations. The black/white statistics is then the 
rate ft = Nf,/N z „. If /3 is plotted versus the coupling parameter p the FigureQ]is obtained. 
The region 0.258 < p < 0.335 where /3 vanishes is remarkable. As stated above, /3 represents 
the rate between the number of black cells and the number of white cells appearing in the 
two-dimensional digitalized representation of the colony evolution. A whole white pattern is 
obtained for this range of p. The phenomenon of spatio-temporal intermittency is displayed 
by the system in the two borders of this parameter region (Fig. [2j. Bursts of low density (black 
color) travel in an irregular way through the high density regions (white color). In this case 
two-dimensional complex patterns are shown by the time evolution of the system (Fig. |2]>c). 
If the coupling p is far enough from this region, i.e., p < 0.25 or p > 0.4, the absorbent region 
loses its influence on the global dynamics and less structured and more random patterns than 
before are obtained (Fig. |2^-d). For p = we have no coupling of the maps, and each map 
generates so called fully developed chaos, where the invariant measure is well-known to be 




Fig. 4. (•) Cjg versus p. The peaks of this modified LMC complexity are also evident just on 
the borders of the absorbent region 0.258 < p < 0.335, where /3 = (x). (Computations have 
been performed with Ap = 0.005 for a lattice of 10000 sites after a transient of 5000 iterations 
and a running of other 2000 iterations). 

symmetric around 0.5. From this we conclude that fi(p = 0) = 1. Let us observe that this 
symmetrical behavior of the invariant measure is broken for small p, and f> decreases slightly 
in the vicinity of p = 0. 

If the LMC complexi ties are quantified as f unction of p, our intuition is confirmed. The 
method proposed in ^Lopez-Ruiz et all Il995h to calculate C is now adapted to the case of 
two-dimensional patterns. First, we let the system evolve until the asymptotic regime is at- 
tained. This transient is discarded. Then, for each time n, we map the whole lattice in a binary 
sequence: if x" < 0.5 and 1 if x" > 0.5, for i = 1, . . . , L. This L-binary string is analyzed by 
blocks of n bits, where n can be considered the scale of observation. For this scale, there are 
2"° possible states but only some of them are accessible. These accessible states as well as their 
probabilities are found in the L-binary string. Next, the magnitudes H, D, Djg, C and Cjg are 
directly calculated for this particular time n by applying the formulas jll3\ . We repeat this 
process for a set of successive time units (n,n + 1, ■ ■ ■ ,w + m). The mean values of H, D, Djg, 
C and Cjg f or these m time units are finally obtained and plotted in Fig. 13141 
Figures l3l4l show the result for the case of n = 10. Let us observe that the highest C and Cjg are 
reached when the dynamics displays spatio-temporal intermittency, that is, the most complex 
patterns are obtained for those values of p that are located on the borders of the absorbent 
region 0.258 < p < 0.335. Thus the plot of C and Cjg versus p shows two tight peaks around 
the values p = 0.256 and p = 0.34 (Fig. 13141 . Let us remark that the LMC complexity C can be 
neglected far from the absorbent region. Contrarily to this behavior, the magnitude Cjg also 
shows high peaks in some other sharp transition of f> located in the region < p < 25, and an 
intriguing correlation with the black/white statistics in the region 0.4 < p < 1. All these facts as 



well as the stability study of the different dynamical regions of system (|4j are not the object of 
the present writing but they deserve attention and a further study. 

If the detection of complexity in the two-dimensional case requires to identify some sharp 
change when comparing different patterns, those regions in the parameter space where an 
abrupt transition happens should be explored in order to obtain the most complex patterns. 
Smoothness seems not to be at the origin of complexity. As wel l as a selected few distinct 
molecules among all the possible are in the basis of life (|McKayll200^ . discreteness and its 
spiky appearance could indicate the way towards complexity. Let us recall that the distribu- 
tions with the highest LMC complexity are just those distribu tions with a spiky-like appear- 
ance dCalbet & Lopez-Ruizll200ll;lAnteneodo & Plastinol,ll996l) . In this line, the striking result 
here exposed confirms the capability of the LMC-like comple xities for signaling a transition t o 
complex behavior when regarding two-dimensional patterns jSanchez & Lopez-Ruizll2005-bh . 

3. Detecting Synchronization in Cellular Automata by Complexity Measurements 

Despite all the efforts devoted to understand the meaning of complexity, we still do not have 
an instrument in the laboratories specially designed for quantifying this property. Maybe this 
is not the final objective of all those theoretical attempts carried out in the most diverse fields 
of knowledge in the last years teolmogorovlll965l:IChaitinl.ll966l:lLempel & Zivl.ll976l:lBenneti 
ll985l;lGrassberger[|l986l;lLloyd & Pagelsill988l;ICruthfieldr989lJshiner et al.l,ll999h , but, for a 
moment, let us think in that possibility. 

Similarly to any other device, our hypothetical apparatus will have an input and an output. 
The input could be the time evolution of some variables of the system. The instrument records 
those signals, analyzes them with a proper program and finally screens the result in the form 
of a complexity measurement. This process is repeated for several values of the parameters con- 
trolling the dynamics of the system. If our interest is focused in the most complex configuration 
of the system we have now the possibility of tuning such an state by regarding the complexity 
plot obtained at the end of this process. 

As a real applicability of this proposal, let us apply it to an a-la-mode problem. The clusteriza- 
tion or sync hronization of chaotic coupled elements was pu t in evidence at the beginning of 
the nineties l|KanekoLll989l;lLopez-Ruiz & Perez-GarciaLll99lh . Since then, a lot of publications 
have been devoted to this subject jBoccaletti et all |2002|) . Let us consider one particular of 
these systems to illuminate our proposal. 

(1) SYSTEM: We take two coupled elementary one dimensional cellular automata (CA: see 
next sec tion in which C A are concisely explained) displaying complex spatio-temporal dy- 
namics dWolframlll983l). It has been shown that this system can undergo through a synchro- 
nization transition jMorelli & Zanettelll998h . The transition to full synchronization occurs at 
a critical value p c of a synchronization parameter p. Briefly the numerical experiment is as 
follows. Two L-cell CA with the same evolution rule <t> are started from different random 
initial conditions for each automaton. Then, at each time step, the dynamics of the coupled 
CA is governed by the successive application of two evolution operators; the independent 
evolution of each CA according to its corresponding rule <3> and the application of a stochastic 
operator that compares the states of and of of all the cells, i = 1,...L, in each automaton. If 
of = of , both states are kept invariant. If of ^ of, they are left unchanged with probability 
1 — p, but bo th states are updated eith er to of or to of with equal probability p II. It is shown 
in reference iMorelli & Zanettel, Il998l) that there exists a critical value of the synchronization 
parameter (p c = 0.193 for the rule 18) above for which full synchronization is achieved. 




Rule 90 Rule 110 




Fig. 5. Spatio-temporal patterns just above the synchronization transition. The left and the 
right plots show 250 successive states of the two coupled automata and the central plot is the 
corresponding difference automaton for the rules 22, 30, 90 and 110. The number of sites is 
L = 100 and the coupling probability is p = 0.23. 



(2) DEVICE: We choose a particular instrument t o perform our measurem ents, that is capable 
of displaying the value of the LMC complexity (C) ILopez-Ruiz et al.l,ll995h defined as in Eq. (TJ, 
C({p;}) = H({Pi}) ' D({pj}) , where {p/} represents the set of probabilities of the N accessible 
discrete states of the system, with pi > , i = 1, • • • , N, and k is a constant. If k = 1/ logN then 
we have the normalized complexity. C is a statistical measure of complexity that identifies 
the entropy or information stored in a system and its disequilibrium, i.e., the distance from 
its actual state to the probability distribution of equilibrium, as the two basic ingredients for 
calculating the complexity of a system. This quantity vanishes both for completely ordered 
and for completely random systems giving then the correct asymptotic properties required 
for a such well-behaved measure, and its calculation has been useful to successfully discern 
many situations regarded as complex. 

(3) INPUT: In particular, the evolution of two coupled CA evolving under the rules 22, 30, 90 
and 110 is analyzed. The pattern of the difference automaton will be the input of our device. 
In Fig. |5]it is shown for a coupling probability p = 0.23, just above the synchronization transi- 
tion. The left and the right plots show 250 successive states of the two automata, whereas the 
central plot displays the corresponding difference automaton. Such automaton is constructed 
by comparing one by one all the sites (L = 100) of both automata and putting zero when the 
states crj and of, i = 1 . . . L, are equal or putting one otherwise. It is worth to observe that 
the difference automaton shows an interesting complex structure close to the synchronization 
transition. This complex pattern is only found in this region of parameter space. When the 
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Fig. 6. The normalized complexity C versus the coupling probability p for different scales 
of observation: n = 1 (°),n = 4 (*) and n = 6 (□). C has been calculated over the last 300 
iterations of a running of 600 of them for a lattice with L = 1000 sites. The synchronization 
transition is clearly depicted around p « 0.2 for the different rules. 



system if fully synchronized the difference automaton is composed by zeros in all the sites, 
while when there is no synchronization at all the structure of the difference automaton is com- 
pletely random. 

(4) METHOD OF MEASUREMENT: How to perform the measure ment of C for such two- 
dimen sional patterns has been presented in the former section jSanchez & Lopez-RuizL 
l2005-ah . We let the system evolve until the asymptotic regime is attained. The variable of 
in each cell of the difference pattern is successively translated to an unique binary sequence 
when the variable i covers the spatial dimension of the lattice, i = l,...,L, and the time vari- 
able 11 is consecutively increased. This binary string is analyzed in blocks of n bits, where n 
can be considered the scale of observation. The accessible states to the system among the 2"° 
possible states is found as well as their probabilities. Then, the magnitudes H, D and C are 
directly calculated and screened by the device. 

(5) OUTPUT: The results of the measurement are shown in Fig. [6] The normalized com- 
plexity C as a function of the synchronization parameter p is plotted for different coupled 
one-dimensional CA that evolve under the rules 22, 30, 90 and 110 , which are known to gen- 
erate complex patterns. All the plots of Fig. [6]were obtained using the following parameters: 
number of cell of the automata, L = 1000; total evolution time, T = 600 steps. For all the cases 
and scales analyzed, the statistical complexity C shows a dramatic increase close to the syn- 
chronization transition. It reflects the complex structure of the difference automaton and the 



capability of the measurement device here proposed for clearly signaling the synchronization 
transition of two coupled CA. 

These results are in agreement with the measurements of C performed in the pat- 
terns generated by a one-dim ensional logistic coupled map lattice in the former section 
jSanchez & Lopez-Ruizt l2005-al) . There the LMC statistical complexity (C) also shows a sin- 
gular behavior close to the two edges of an absorbent region where the lattice displays spatio- 
temporal intermittency. Hence, in our present case, the synchronization region of the coupled 
systems can be interpreted as an absorbent region of the difference system. In fact, the highest 
complexity is reached on the border of this region for p « 0.2. The parallelism between both 
systems is therefore complete. 

4. Self-Synchronization of Cellular Automata 

Cellular automata (CA) are discrete dynamical systems, discrete both in space and time. The 
simplest one dimensional version of a cellular automaton is formed by a lattice of N sites or 
cells, numbered by an index i = 1,...,N, and with periodic boundary conditions. In each site, 
a local variable c, taking a binary value, either or 1, is asigned. The binary string c(f) formed 
by all sites values at time t represents a configuration of the system. The system evolves in 
time by the application of a rule O. A new configuration a(t + 1) is obtained under the action 
of the rule O on the state c(f). Then, the evolution of the automata can be writen as 

<r(f + l) = *[(r(f)]. (6) 

If coupling among nearest neighbors is used, the value of the site i, tr,-(f + 1), at time t + 1 is 
a function of the value of the site itself at time t, cr,-(f), and the values of its neighbors <7y_i(f) 
and Cj +1 (f) at the same time. Then, the local evolution is expressed as 

o- i (t + l)=<p(cr i - 1 (t),cr i (t),cri +1 (t)), (7) 

being <p a particular realization of the rule <t>. For such particular implementation, there will 
be 2 3 different local input configurations for each site and, for each one of them, a binary value 
can be assigned as output. Therefore there will be 2 8 different rules <p , also called the 'Wolfram 
rules. Each one of these rules produces a different dynamical evolution. In fact, dynamical 
behavior generated by all 256 rules were already classified in four generic classes. The reader 
interested in the details of such classification is addressed to the original reference jWolframl 
Il983h . 

CA provide us with simple dynamical systems, in which we would like to essay differ- 
ent methods_of_s_ynchronization. A stochastic synchronization technique was introduced 
in jMorelli & Zanettel Il998f) that works in synchronizing two CA evolving under the same 
rule "J>. The two CA are started from different initial conditions and they are supposed to have 
partial knowledge about each other. In particular, the CA configurations, u 1 (f ) and o~ 2 (t), are 
compared at each time step. Then, a fraction p of the total different sites are made equal (syn- 
chronized). The synchronization is stochastic since the location of the sites that are going to be 
equal is decided at random. Hence, the dynamics of the two coupled CA, c(f ) = {a 1 (f ), a 2 (f ) ), 
is driven by the successive application of two operators: 

1. the deterministic operator given by the CA evolution rule <t>, <j>[c(f)] = 
(4>[t7 1 (t)],4>[t7- 2 (J)]),and 
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Fig. 7. Rule 90 has two stable patterns: one repeats the Oil string and the other one the 
00 string. Such patterns are reached by the first self-synchronization method but there is a 
dynamical competition between them. In this case p = 0.9. Binary value is represented in 
white and 1 in black. Time goes from top to bottom. 



2. the stochastic operator Tp that produces the result r p [c(f)], in such way that, if the 
sites are different (of ^ of), then Tp sets both sites equal to of with the probability 
p/2 or equal to of with the same probability p/2. In any other case T p leaves the sites 
unchanged. 

Therefore the temporal evolution of the system can be written as 

a(t + l) = (T p0 ®)[a(t)]=r p [(®[a 1 (t)],®[a 2 (t)])]. (8) 

A simple way to visualize the transition to synchrony can be done by displaying the evolution 
of the difference automaton (DA), 

S,(t)=\o}(t)-of(t)\. (9) 
The mean density of active sites for the DA 

i N 
(=1 

represents the Hamming distance between the automata and verifies < p < 1. The automata 
will b e synchronized when limt_ s . 00| o(f) = 0. As it has been described in jMorelli & Zanettel 
Il998l) that two different dynamical regimes, controlled by the parameter p, can be found in 
the system behavior: 

p < p c —5- limf-yoo p(0 7^ (no synchronization), 
p > p c —5- limf_ ! , 00 p(i) = (synchronization), 
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Fig. 8. Mean density p vs. pmax = p for different rules evolving under the second synchro- 
nization method. The existence of a transition to a synchronized state can be clearly observed 
for rule 18. 



being p c the parameter for which the transition to the synchrony occurs. When p < p c complex 
structures can be observed in the DA time evolution. In Fig. [5] typical cases of such behavior 
are shown near the synchronization transition. Lateral panels represent both CA evolving in 
time where the central strip displays the evolution of the corresponding DA. When p comes 
close to the critical value p c the evolution of S(t) becom es rare and resembles the problem of 
structures trying to percolate in the plane jPomeaul.ll986|) . A method to detect this kind of tran- 
sition, based in the calculatio n of a statistical measure of com p lexity for patterns, has been pro - 
posed in the former sections dSanchez & Lopez-Ruizll2005-al) . l lSanchez & Lopez-Ruizll2005-bh . 

4.1 First Self-Synchronization Method 

Let us now take a single cellular automaton (ToffoU & Margolusl Il987l; lllachinskil EzOOlh . If 
c 1 (f) is the state of the automaton at time t, a 1 (f) = c(f ), and c 2 (f) is the state obtained from 
the application of the rule O on that state, o 2 (f ) = Ofu 1 (f )], then the operator T p can be applied 
on the pair (a l {t),a 2 (t)), giving rise to the evolution law 



updated to the correspondent values taken in cr 1 (i) with a probability p. The updated array 
c 2 (t) is the new state a(t + 1). 

It is worth to observe that if the system is initialized with a configuration constant in time 
for the rule <E>, 4>[<r] = a, then this state a is not modified when the dynamic equation fiTi is 
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Fig. 9. Mean density p vs. p for rule 18 evolving under the third self-synchronization method. 
The existence of a transition to a synchronized state can be observed despite of the random- 
ness in the election of neighbors within a range L, up to L = 4. 

applied. Hence the evolution will produce a pattern constant in time. However, in general, 
this stability is marginal. A small modification of the initial condition gives rise to patterns 
variable in time. In fact, as the parameter p increases, a competition among the different 
marginally stable structures takes place. The dynamics drives the system to stay close to 
those states, although oscillating continuously and randomly among them. Hence, a complex 
spatio-temporal behavior is obtained. Some of these patterns can be seen in Fig. [7] However, 
in rule 18, the pattern becomes stable and, independently of th e initial conditions, the system 
evolves toward this state, which is the null pattern in this case dSanchez & Lopez-Ruizll2006l) . 

4.2 Second Self-Synchronization Method 

Now we introduce a new stochastic element in the application of the operator Tp. To differ- 
entiate from the previous case we call it Tp. The action of this operator consists in applying at 
each time the operator Tp, with p chosen at random in the interval (0, p). The evolution law 
of the automaton is in this case: 

tr(f + 1) = TpKo-HtWit))} = Tpl(<r(t)M<r(t)})}- (12) 

The DA density between the present state and the previous one, defined as S(t) = \ c(f) — a(t — 
1) , is plotted as a function of p for different rules $ in Fig. [8] Only when the system becomes 
self-synchronized there will be a fall to zero in the DA density. Let us observe again that the 
behavior reported in the first self-synchronization method is newly obtained in this case. Rule 
18 undergoes a phase transition for a critical value of p. For p greater than the critical value, 




Fig. 10. Mean density p vs. p for different rules evolving under the third self-synchronization 
method. The density of the system decreases linearly with p. 



the method is able to find the stable structure of the system iSanchez & Lopez-RuizLl2006h . For 
the rest of the rules the freezing phase is not found. The dynamics generates patterns where 
the different marginally stable structures randomly compete. Hence the DA density decays 
linearly with p (see Fig. |8). 

4.3 Third Self-Synchronization Method 

At last, we introduce another type of stochastic element in the application of the rule 4>. Given 
an integer number L, the surrounding of site i at each time step is redefined. A site i; is 
randomly chosen among the L neighbors of site i to the left, (i — L,...,i — 1). Analogously, a 
site zV is randomly chosen among the L neighbors of site i to the right, (z + 1, ...,/ + L). The rule 
<I> is now applied on the site i using the triplet (i;,i,z r ) instead of the usual nearest neighbors 
of the site. This new version of the rule is called being <J>l=i = ®- Later the operator Tp 
acts in identical way as in the first method. Therefore, the dynamical evolution law is: 

cr(f + 1) = T v [{a\t),a\t))\ = T p [(<r(t),<P L [cr(t)})}. (13) 

The DA density as a function of p is plotted in Fig. [9] for the rule 18 and in Fig. [l0]for other 
rules. It can be observed again that the rule 18 is a singular case that, even for different L, 
maintains the memory and continues to self -synchronize. It means that the influence of the 
rule is even more important than the randomness in the election of the surrounding sites. The 
system self-synchronizes and decays to the corresponding stable structure. Contrary, for the 
rest of the rules, the DA density decreases linearly with p even for L = 1 as shown in Fig. [lOl 
The systems oscill ate randomly among their di fferent marginally stable structures as in the 
previous methods {Sanchez & Lopez-Ruizl.l2006|) . 
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Fig. 11. Space-time configurations of automata with N = 100 sites iterated during T = 400 
time steps evolving under rules 18 and 150 for p < p c . Left panels show the automaton evolu- 
tion in time (increasing from top to bottom) and the right panels display the evolution of the 
corresponding DA. 

5. Symmetry Pattern Transition in Cellular Automata with Complex Behavior 

In this section, the stoc hastic synchronization method introduced in the former sections 
jMorelli & Zanettd. Il998h for two CA is specifically used to find symmetrical patterns in the 
evolution of a single automaton. To achieve this goal the stochastic operator, below described, 
is applied to sites symmetrically located from the center of the lattice. It is shown that a sym- 
metry transition take place in the spatio-temporal pattern. The transition forces the automaton 
to evolve toward complex patterns that have mirror symmetry respect to the central axe of 
the pattern. In consequence, this synchronization method can also be interpreted as a control 
technique for stabilizing complex symmetrical patterns. 

Cellular automata are extended systems, in our case one-dimensional strings composed of N 
sites or cells. Each site is labeled by an index i = 1, . . . , N, with a local variable s,- carrying a 
binary value, either or 1. The set of sites values at time t represents a configuration (state 
or pattern) at of the automaton. During the automaton evolution, a new configuration c t+1 at 
time t + 1 is obtained by the application of a rule or operator 4> to the present configuration 
(see former section): 

a t+l = O [a t \ . (14) 



Rule 18 



Rule 150 




Fig. 12. Time configurations of automata with N = 100 sites iterated during T = 400 time steps 
evolving under rules 18 and 150 using p > p c . The space symmetry of the evolving patterns is 
clearly visible. 
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Fig. 13. Asymptotic density of the DA for different rules is plotted as a function of the coupling 
probability p. Different values of p c for each rule appear clearly at the points where p -> 0. 
The automata with N = 4000 sites were iterated during T = 500 time steps. The mean values 
of the last 100 steps were used for density calculations. 

5.1 Self-Synchronizati on Method by Symmetr y 

Our present interest ISanchez & Lopez-Ruia. I2008D resides in those CA evolving under 
rules capable to show asymptotic complex behavior (rules of class III and IV). The tech- 
nique a pplied here is similar to the synchronization scheme introduced by Morelli and 
Zanette {Morelli & Zanettelll998h for two CA evolving under the same rule <E>. The strategy 
supposes that the two systems have a partial knowledge one about each the other. At each time 
step and after the application of the rule <S>, both systems compare their present configurations 
<&[c}] and <!>[<T f 2 ] along all their extension and they synchronize a percentage p of the total of 
their different sites. The location of the percentage p of sites that are going to be put equal is 
decided at random and, for this reason, it is said to be an stochastic synchronization. If we call 
this stochastic operator Tp, its action over the couple (3>[(Tj],3>[(7^]) can be represented by the 
expression: 

(oi +1 ,of +1 )=T p (<l>[o}],9[°t\) = {T p ° of). (15) 

The same strategy can be a pplied to a single automaton with a even number of sites 
JSanchez & Lopez-R uiz, 2008). Now the evolution equation, <r f+1 = (T p o <J>)[cf], given by 
the successive action of the two operators O and Tp, can be applied to the configuration Of as 
follows: 

1. the deterministic operator $ for the evolution of the automaton produces <E>[<7}], and, 

2. the stochastic operator Tp, produces the result T p (<& [cr t ] ), in such way that, if sites sym- 
metrically located from the center are different, i.e. s, ^ s N _ !+ i, then T p equals s N _,- +1 
to Sj with probability p. Tp leaves the sites unchanged with probability 1 — p. 
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54 


60 
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110 
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146 


150 


182 
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0.25 


0.27 


1.00 


0.20 


1.00 


0.25 


0.37 


1.00 


0.27 


0.30 


0.25 


0.37 


0.25 



Table 1. Numerically obtained values of the critical probability p c for different rules displaying 
complex behavior. Rules that can not sustain symmetric patterns need fully coupling of the 
symmetric sites, i.e. (p c = 1). 



A simple way to visualize the transition to a symmetric pattern can be done by splitting the 
automaton in two subsystems (of, of), 

• o} , composed by the set of sites s (i ) with / = 1, . . . , N/2 and 

• of, composed the set of symmetrically located sites s(N — i + 1) with / = 1,.. .,N/2, 
and displaying the evolution of the difference automaton (DA), defined as 

S'=\o}-of\. (16) 

The mean density of active sites for the difference automaton, defined as 

o N/2 

1=1 

represents the Hamming distance between the sets o l and o 2 . It is clear that the automaton 
will display a symmetric pattern when Iirrtf-^oop* = 0. For class III and IV rules, a symmetry 
transition controlled by the parameter p is found. The transition is characterized by the DA 
behavior: 

when p < p c — > \m\ t ^cop t (complex non-symmetric patterns), 
when p > p c — > lirrv^co/c/ = (complex symmetric patterns). 

The critical value of the parameter p c signals the transition point. 

In Fig. [ll] the space-time configurations of automata evolving under rules 18 and 150 are 
shown for p < p c . The automata are composed by N = 100 sites and were iterated during 
T = 400 time steps. Left panels show the automaton evolution in time (increasing from top 
to bottom) and the right panels display the evolution of the corresponding DA. For p < p c , 
complex structures can be observed in the evolution of the DA. As p approaches its critical 
value p c , the evolution of the DA become more s t umped and reminds the proble m of struc- 
tures trying to percolate the plane fPomeau, 19861: ISanchez & Lopez-Ruizll2005-a|) . In Fig. [12] 
the space-time configurations of the same automata are displayed for p > p c . Now, the space 
symmetry of the evolving patterns is clearly visible. 

Table[l]shows the numerically obtained values of p c for different rules displaying complex be- 
havior. It can be seen that some rules can not sustain symmetric patterns unless those patterns 
are forced to it by fully coupling the totality of the symmetric sites (p c = 1). The rules whose 
local dynamics verify (p (si,sq,S2) = <P(s2,sq,s\) can evidently sustain symmetric patterns, and 
these structures are induced for p c < 1 by the method here explained. 

Finally, in Fig. [I3]the asymptotic density of the DA, p t for t — > oo, for different rules is plotted 
as a function of the coupling probability p. The values of p c for the different rules appear 
clearly at the points where p — > 0. 



6. Conclusion 



A method to measure statistical complexity in extended systems has been implemented. It 
has been applied to a transition to spatio-temporal complexity in a coupled map lattice and 
to a transition to synchronization in two stochastically coupled cellular automata (CA). The 
statistical indicator shows a peak just in the transition region, marking clearly the change of 
dynamical behavior in the extended system. 

Inspired in stochastic synchronization methods for CA, different schemes for self- 
synchronization of a single automaton have also been proposed and analyzed. Self- 
synchronization of a single automaton can be interpreted as a strategy for searching and con- 
trolling the structures of the system that are constant in time. In general, it has been found 
that a competition among all such structures is established, and the system ends up oscillat- 
ing randomly among them. However, rule 18 is a unique position among all rules because, 
even with random election of the neighbors sites, the automaton is able to reach the configu- 
ration constant in time. 

Also a transition from asymmetric to symmetric patterns in time-dependent extended systems 
has been described. It has been shown that one dimensional cellular automata, started from 
fully random initial conditions, can be forced to evolve into complex symmetrical patterns by 
stochastically coupling a proportion p of pairs of sites located at equal distance from the center 
of the lattice. A nontrivial critical value of p must be surpassed in order to obtain symmetrical 
patterns during the evolution. This strategy could be used as an alternative to classify the 
cellular automata rules -with complex behavior- between those that support time-dependent 
symmetric patterns and those which do not support such kind of patterns. 
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